Finite DifFerence-Time Domain solution of Dirac equation and the 

Klein Paradox 



Neven Simicevic 

Center for Applied Physics Studies, Louisiana Tech University, Ruston, Louisiana 71272, USA 



Abstract 

The time-dependent Dirac equation is solved using the three-dimensional Finite Difference-Time 
Domain (FDTD) method. The dynamics of the electron wave packet in a scalar potential is 
studied in the arrangements associated with the Klein paradox: potential step barriers and linear 
potentials. No Klein paradox is observed. 
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In our previous paper Q, the Finite Diff'erence-Time Domain (FDTD) method, originally 
introduced by Kane Yee (T\ to solve Maxwell's equations, was for the first time applied to solve 
the three-dimensional Dirac equation. The Zitterbewegung and the dynamics of a well-localized 
electron were used as examples of FDTD applied to the case of free electrons. In this paper the 
Finite Diff'erence-Time Domain method is applied to the case of the motion of an electron wave 
packet inside and scattering from the potential step barrier or linearly dependent potential. 

In the analysis of the Dirac equation, Oskar Klein calculated that if the electrons scatter from 
the repulsive potential step of sufficient strength, they were partially and without impedance 
transmitted even if the transmission was forbidden by the conservation of energy [3 |. At the 
same time, it was possible that the reflected current was larger than the incoming one. Since the 
consequences of such a finding are counterintuitive and lead to interpretational problems of the 
Dirac equation, the effect is known as the Klein paradox. For comparison, if the same scattering 
was analyzed in the non-relativistic limit by solving the Schrodinger equation, the total reflection 
of electrons was obtained, as intuitively expected. 

Despite its counterintuivity, the Klein paradox profoundly influenced the development of 
relativistic quantum mechanics becoming one of its fundamentals. It facilitated the transition 
from the single-particle interpretation of the Dirac equation into interpreting the Dirac field as 
a many -body problem |4|. Much work in nuclear, particle, and astro-physics, not all consistent 
with each other, has been published on the Klein paradox. A recent historical review of the 
paradox was written by Dombey and Calogeracos |5 |. Krekora et al, in addition to publishing 
the results on their temporally resolved numerical solution of the Dirac equation, provided the 
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status of recent theoretical investigation |[6j|7l|8j|9l. While most of the work conforms with the 
Klein solution and the interpretation of the Klein paradox, there is still some dissent, with the 
most recent contained in Ref . (TOl [TT] [121 • Recently, the interest in the Klein paradox increased 
with the proposal to experimentally test it in a simple experiment using electrostatic barriers in 
single- and bi-layer graphene (TSj, but until now, no conclusive experimental evidence of the 
Klein paradox has been found ifTTIl . 

The Klein paradox results from the solution of the stationary time-independent Dirac equa- 
tion. In this work, the results of the solution of the time-dependent Dirac equation are presented. 
The FDTD scheme is applied to the Dirac equation for the case when the electromagnetic field 
described by the four-potential = {Aq(x),A(x)} is minimally coupled to the particle 141 [T4l 



ih 



~dt 



where 



and 



Hfree = -icHa • V -h /3mc 
Hint = -ea • A -h Mo, 



( ^iW ^ 



(1) 

(2) 
(3) 

(4) 



The matrices a and jS are expressed using 2x2 Pauli matrices cr's and the 2 x 2 unit matrix /. 

The FDTD schematics to solve Eq. ([T]) follows Yee's leapfrog algorithm |1, 2|. The wave 
functions and ^2 at the time ^ - 1/2 are used to calculate the wave functions ^3 and ^4 at 
the time n, which are then used to calculate the wave functions ^1 and ^2 at the time n -\- 1/2, 
and so on. The same numerical requirements as in the case of electrodynamics are followed 111. 
While the dynamics of a Dirac electron can be studied in an environment described by any four- 
potential A^ regardless of its complexity and time dependency, the study of the Klein paradox 
requires only that Ao(x) 9^ in the designated region and A(x) = everywhere. 

As a consequence of the Dirac equation being of first order and linear in d/dt, the entire 
dynamics of the electron is defined, as in the case of Maxwell's equations, only by its initial 
wave function. The dynamics of a wave packet is defined by its initial wave function 
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where N = [(27t)^^^Xq]~^^^ is normalizing constant. Eq. ([sj) represents a wave packet whose initial 
probability distribution is of a normalized Gaussian shape. Its size is defined by the constant xq, 
its spin is pointed along the z-axis, and its motion is defined by the values of pi,p2, and 773. In 
the "single particle interpretation" of the Dirac equation, if we choose P2 = P3 = and pi 
the wave packet should move in the x-direction. This is not the case. Because of the localization 
of the wave packet and limitation of the direction of the spin, the initial condition in Eq. ^ may 
contain also a component moving in the opposite direction flJIISl- For the case of free electrons 
some of the dynamics of this wave packet were studied in Ref. Q. 
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Figure 1: Six stages of motion of the wave packet initialized by Eq. (|5} in the x-z plane. After splitting in two as a result 
of the initial conditions, the larger component of the wave packet moves in the +x-direction and completely reflects from 
the repulsive supercritical potential step barrier. One can observe the interference between incoming and reflecting wave, 
but no Klein paradox. The animation can be accessed on-line fTTI . 



The dynamics of the wave packet is very complex. The energy of the particle described by 
the packet depends on its localization, defined by the Gaussian component of the wave function, 
and its initial momentum, part of the wave function's phase. While an extensive study was done 
in order to understand the interplay between the height and the shape of the potential barrier, and 
particle localization and its initial momentum, in this paper, we report only the results of several 
cases associated with the Klein paradox: the penetration into a potential barrier of supercritical 
potential satisfying the condition eV > E + mc^. 

The first part of this paper describes the dynamics of scattering of the wave packet initiated 
by Eq. ^ from the potential barrier 

V for jc > 
forx<0 ' 

where x represents the coordinate in the x-direction. The initial values of the wave packet mo- 
menta were pi = 18.75M^V/c and P2 = P3 = 0, and the size was defined by xq = 10"^^ m. The 
step barrier repulsive potential was V = 25 x 10^ V, satisfying eV > E -\- mc^. The dynamics of 
the scattering of this wave packet from the repulsive potential is shown in Fig. [T] After splitting 
in two as a result of the initial conditions |[I|[l5][T6J, the larger component of the wave packet 
moves in the direction of the repulsive supercritical potential step barrier and completely reflects 
from it. While we can observe the interference between the incoming and reflected waves in 
front of the potential barrier, no Klein paradox is observed. The wave packet did not penetrate 
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Figure 2: Six stages of motion of the wave packet initialized by Eq. |5} in the x-z plane. The wave packet partially 
reflects and partially penetrates the attractive supercritical potential step barrier. In addition to the interference between 
incoming and reflecting wave, the shape of the wave packet changes as it penetrates the barrier. The animation can be 
accessed on-line fTTl . 

the potential barrier as predicted by the Klein solution of the Dirac equation, but completely re- 
flected back as in the case of non-relativistic quantum mechanics. For comparison, the dynamics 
associated with the attractive potential step barrier of the same magnitude is shown in Fig. [2] 
The wave packet partially reflects and partially penetrates the potential step barrier. Integration 
of the probability density function shown in Fig. [s] gives the reflection and transmission 
coefficients R and T. In the case of the repulsive supercritical potential step barrier, as already 
discused, R = I and T = 0. (Increasing the step barrier repulsive potential did not change this 
result.) In the case of the attractive supercritical potential step barrier R = 0.46 and T = 0.54. 

The cases of attractive or repulsive subcritical potential barriers, when eV < E -\- mc^, also 
show nonpeculiar behavior. The same wave packet was scattered from the step barrier potential 
of y = ±50 X 10^ V. The dynamics of the scattering from the repulsive potential is shown in 
Fig. [4] and from the attractive potential in Fig. [5] In the case of the repulsive potential after the 
interference at the potential barrier, part of the wave packet is reflected and part penetrates the 
barrier. The wave packet which has penetrated the potential barrier broadened and dissipated, 
showing, as expected in the non-relativistic case, damping with distance. In the case of the at- 
tractive potential part of the wave packet which has penetrated the potential barrier showed little 
broadening, comparable to the distortion of the free electron wave packet 1 1 1. The shapes of the 
wave packets after the penetration into the potential barriers are shown in Fig. [6] To test the 
efl'ect of the shape of the potential barrier and the efl'ect of the wave packet shape and momentum, 
two more tests were performed. To verify that the properties of the wave packet propagation are 
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Figure 3: The probability density function for the dynamics associated with the repulsive supercritical potential step 
barrier (left) and with the attractive supercritical potential step barrier (right). is shown at the time before the wave 
packet interacts with potential barrier (top), at the time of the interaction (in the middle), and at the time when the wave 
packet is completely reflected from the repulsive barrier or is completely separated for the attractive barrier (bottom). 




Figure 4: Six stages of motion of the wave packet initialized by Eq. j5j in the x-z plane for the case of the repulsive 
subcritical potential step barrier. After the interference stage, the major portion of the wave packet is reflected while a 
small part penetrates the barrier and dissipates. The animation can be accessed on-line 1 17 1. 
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Figure 5: Three stages of motion of the wave packet initialized by Eq. j5} in the x-z plane for the case of the attractive 
subcritical potential step barrier. This time, after the interference stage, a major portion penetrates the barrier. The 
animation can be accessed on-line UTi . 




Figure 6: Four time steps of the probability density function for the dynamics associated with the repulsive (left) and 
the attractive (right) subcritical potential step barrier. In the case of the repulsive potential, broadens and dissipates 
as it penetrates the potential barrier. In the case of the attractive potential, the distortion of is small, comparable to 
the distortion of a free electron wave packet [IJ. 
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Figure 7: Three stages of motion of the wave packet initialized by Eq. |5} in the x-z plane for the case of the repulsive 
potential barrier described by Eq. Left and right figures correspond to the wave packet before and after the scattering. 
Middle figure corresponds to the maximal penetration. The lines correspond to the region a > x > in Eq. JTJ. The 
animation can be accessed on-line i ifTI . 
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Figure 8: Three stages of motion, in the x-z plane, of the wave packet of the momentum pi = 0.1S15MeV/c and 
P2 = P3 = ^ initialized by Eq. js}. Top row is for the supercritical and bottom row for the subcritical repulsive potential. 
Left and right figures correspond to the wave packet before and after the scattering. Middle figure corresponds to the 
time of interaction of the wave packet with the potential barrier. The line corresponds to x = 0. The animation can be 
accessed on-line L17J . 
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Figure 9: Three stages of motion of the wave packet in the x-z plane for the case of the angle of incidence of 45°. Top 
row is for the repulsive and the bottom row for attractive supercritical potential step barrier. In the case of the repulsive 
potential the wave packet reflects from the potential with the reflected angle equal to the incident angle. In the case of 
the attractive potential, part of the wave packet reflects from and part refracts into the potential. The animation can be 
accessed on-line 1 17.1 . 



not due to the shape of the potential barrier, the step potential barrier, Eq. ([6]), was replaced by 
the potential representing a constant electric field in a finite region of space 

IV for X > a 
-EX for a > X > . (7) 
for jc < 

In this particular case, a = 50 x 10"^"^ m, e = -50 x 10^^ V/m, and V = 25 x 10^ V. x represents 
the coordinate in the x-direction. The snapshots of the dynamics of the scattering of the same 
wave packet, Eq. Q, from the repulsive potential of this form are shown in Fig. [7] Again, the 
wave packet did not penetrate the potential and no Klein paradox was observed. 

The importance of the wave packet momentum and shape was tested by reducing the initial 
momentum to p\ - 0.1875M^V/c and pi - - 0, and increasing the width of the wave 
packet to = 10"^^m. This reduced the momentum due to the particle localization to A/? ~ 
fi/Ax ~ 0.2M^V/c, keeping the energy of the wave packet below particle production threshold 
and preventing generation of appreciable amounts of negative-energy components of the wave 
packet |1, 14 1 . The step barrier supercritical repulsive potential, satisfying eV > E ^- mc^, was 
chosen V = 1.5 x 10^ V, and subcritical repulsive potential, satisfying eV < E-\-mc^, was chosen 
V = 1.0 X 10^ V. The snapshots of the dynamics of the scattering in both cases are shown in 
Fig. [8] The wave packet did not penetrate the supercritical repulsive potential and no Klein 
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paradox was observed. The wave packet did penetrate the subcritical repulsive potential with the 
reflection and transmission coefficients ofR^ 0.6 and T ^ 0.4 

We also tested the eff'ect of the angle of impact. The initial values of the wave packet mo- 
menta were Pi = P3 = l^.l5MeV/c and p2 = 0, and the step barrier repulsive potential was 
y = 25 X 10^ V. Such a wave packet scatters from the potential at an angle of 45°. The dynamics 
of this scattering is shown in Fig. |9] The wave packet reflects from the potential similarly to 
the reflection of the electromagnetic wave with the reflected angle equal to the incident angle. 
The wave packet did not penetrate the potential barrier and, again, the Klein paradox was not ob- 
served. When the repulsive potential was replaced with an attractive potential, part of the wave 
packet reflected and part refracted into the potential. This scattering is also shown in Fig. [9] 

Finally, to test spin orientation, some of the studies were repeated for the dynamics of the 
spin-flipped wave packet defined by its initial wave function 



The results did not change. 

To conclude, the full three-dimensional Finite Diff'erence Time Domain (FDTD) method was 
developed to solve the Dirac equation. In this paper, the method was applied to the dynamics 
of a Dirac electron in a scalar potential, particularly to the arrangements corresponding to the 
dynamics associated with the Klein paradox. The Klein paradox, resulting from the solution of 
the stationary Dirac equation, predicts unimpeded penetration of particles in the energy forbid- 
den region, in contradiction with intuition and the behavior described by Schrodinger equation. 
Solutions of the time-dependent Dirac equation presented in this paper show no such penetration 
leading to the conclusion that the Klein paradox does not exist. 
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acknowledged. 
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